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Abstract. Polynomial chaos expansions (PCE) are widely used in the framework of uncertainty 
quantification. However, when dealing with high dimensional complex problems, challenging 
issues need to be faced. For instance, high-order polynomials may be required, which leads to a 
large polynomial basis whereas usually only a few of the basis functions are in fact significant. 
Taking into account the sparse structure of the model, advanced techniques such as sparse 
PCE (SPCE), have been recently proposed to alleviate the computational issue. In this paper, 
we propose a novel approach to SPCE, which allows one to exploit the model’s hierarchical 
structure. The proposed approach is based on the adaptive enrichment of the polynomial basis 
using the so-called principle of heredity. As a result, one can reduce the computational burden 
related to a large pre-defined candidate set while obtaining higher accuracy with the same 
computational budget. 
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1 INTRODUCTION 

Nowadays, uncertainty quantification (UQ) has become a key topic in computational science 
and engineering applications. Its formulation comprises uneertainty propagation, whieh deals 
with propagating the uncertainties in the input parameters of a computational model to the 
output quantities of interest. Typical applications include, e.g. eonducting reliability analysis 
of complex systems, estimating the statistical distribution of the outputs, computing sensitivity 
indices, etc. 

Uncertainty propagation is traditionally eonducted by means of sampling-based methods 
sueh as Monte Carlo simulation, importanee sampling and subset simulation. These methods, 
however, are often not suitable for practical problems because of the eomputational eost re¬ 
lated to evaluating complex models repeatedly to obtain a sufficient accuracy. For this reason, 
inexpensive-to-evaluate surrogate models have been of particular interest as substitutes of the 
full models. One way of building surrogate models is by means of the spectral approach, in 
which the quantity of interest (Qol) is eonsidered as a funetion of the input parameters in a 
suitable funetional space and then projeeted onto an appropriate basis. Multiple basis ehoiees 
are available in the literature, including trigonometric, radial basis, wavelet and polynomial 
functions. In particular, the use of polynomial chaos expansions (PCE) has been found highly 
effective in numerous applications, see e.g. llll|2l|3|4||. 

Building a polynomial chaos expansion requires determining the eoefficients corresponding 
to the polynomial ehaos basis. For blaek-box type problems, it is impossible to direetly estimate 
the polynomial chaos coefficients by means of intrusive approaches that require prior knowl¬ 
edge of the governing equations of the systems. In such cases, one has to rely on non-intrusive 
methods such as projection, regression which require only a finite number of observations of 
the system. Among those, least-square minimization methods have received particular interest 
Bil. However, they are faeing the so-called curse-of-dimensionality, i.e. the computational 
cost increases exponentially with the dimension of the input space. This is the ease, e.g. when 
a large number of basis elements is involved but only a small number of experiments (observa¬ 
tions) is available, which is not sufficient for the accurate estimation of the coefficients. Several 
numerieal methods have been proposed recently for alleviating such issue, including compres¬ 
sive sensing Q, / 2 -norm regularized regression IlSl and least angle regression |l9]]. Most of the 
existing methods rely on considering a set of candidate polynomials that is selected a priori, 
then computing the corresponding coefficients either by setting up a constrained optimization 
problem or by making assumptions on the structure of the system. Recently, ifTOll and ifTTll 
proved that in some eases it ean be more effective to apply an adaptive basis PCE approach 
instead. In this paper, we introduee a novel hierarchical method that eonsists in updating the 
candidate set adaptively instead of fixing it. 

The paper is organized as follows: in Section 2, we recall briefly the non-intrusive regression- 
based approach for computing PCE. In Section 3, the hierarchical adaptive method for PCE is 
introduced. Eor the sake of illustration, we prove the effectiveness of the proposed method on 
two numerical examples. 
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2 POLYNOMIAL CHAOS EXPANSIONS 
2.1 Polynomial chaos expansions 

Consider a system with uncertain input parameters whose behavior can be described by the 
following equation: 


Y = M{X), (1) 

in which X = {Xi, ... , Xm}^ is a vector of M independent uncertain input parameters and Y 
is the output quantity of interest. For the sake of simplicity, and without loss of generality, one 
considers only the case of a scalar-valued output. Note that one does not know Eq. ([T]) explicitly 
when the computational model is of a black-box type. 

In the Hilbert space of square-integrable functions of variables {Xi,i = 1, ... , M}^, one 
can select a basis of orthonormal univariate polynomials {V'l, k G N}, also known as polyno¬ 
mial chaos functions ifT^fT^ . associated with the probability density measure (da:*) = /x* (a;*)da;*. 
For instance when X* is a uniform (resp. standard normal) random variable, the corresponding 
basis consists of orthonormal Legendre (resp. Hermite) polynomials. The polynomial chaos 
expansion of Y reads: 


Y = Y, Voci’cciX) (2) 

where a = (oi, ... , ctm) is a multi-index with a*, z = 1, ... , M denoting the degree of the 

M 

polynomial in the direction of X*, ipa{X) = multivariate polynomials obtained 

i 

by the tensor product of univariate polynomials and Ha are the associated coefficients. 

The regression-based adaptive sparse polynomial chaos expansions (SPCE) consists of two 
sequential steps. First, a set of candidate polynomial chaos basis is chosen a priori (see Sec¬ 
tion |2^. Second, the relevant basis are selected from the candidate set and the corresponding 
coefficients are computed (see Section |2i^ . The two steps are described subsequently. 

2.2 Truncation schemes 

In practice, an approximate form of PCE in Eq. Q with a finite number of terms must be 
used: 


Y^Yy<-'^o.{x) ( 3 ) 

ae.4 

The set of multi-indices A is obtained by truncating the space. We present herein two 
truncation schemes that were recently proposed in ifTdll . 

In the first scheme, one assumes that the main effects and the low-degree interaction effects 
are more important than high-degree interaction effects. Thus, in the multi-indices space, the 
relevant terms lie in the subspace defined as follows ifTHl : 

/ M \ 

.A = = {a G : llallg = f j (4) 

where 0 < g < 1 is the parameter governing the hyperbolic truncation scheme and p is the 
prescribed maximum degree of the polynomials. Note that in the special case when q is assigned 
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the value of 1, the standard truncation scheme is recovered: 

M 

A = = {g € : ||q:||i = aj ^ p} (5) 

1=1 

Given a prescribed value of p, the size of the polynomial basis obtained with the hyperbolic 
truncation scheme increases as q increases and attains its maxima at g = 1. 

In the second scheme, one assumes that low-rank effects are more important than high- 
rank effects, i.e. two-dimensional interaction terms are more relevant than three-dimensional 
interaction terms and so on. This assumption is based on the so-called hierarchy principle ifTSll 
which states that the model can be approximated by low-rank terms. The low-rank truncation 
scheme reads El: 

M M 

A = A^’^ = {a e : ||a||o = J]] lai>o < r, ||a||i = ^ p} (6) 

i=l i=l 

in which ||a||o is the rank of the multivariate polynomial 'ipa, defined as the total number of 
non-zero component indices ai,i = 1, ... , M. The prescribed rank r is chosen as a small 
integer value, e.g. r = 2,3, 4. 

The proposed schemes, namely hyperbolic and low-rank truncation, usually result in a can¬ 
didate basis of manageable size. The two strategies consist in selecting a set of candidate basis 
a priori, from which a sparse subset is then retained. The accuracy of the PCE is measured by 
means of error estimates, e.g. empirical error or leave-one-out (LOO) error e^oo El . 

2.3 Least angle regression 

Given a candidate set generated by a truncation scheme (see Section |T2| ) and an experimental 
design (a set of input values X and the corresponding observed output 3 ^ of the model), one has 
to determine the relevant basis and the associated coefficients of the expansion: 

(7) 

o.eA 

This is an example of the variable selection problem in statistics. To this purpose, a wide variety 
of methods is available, including ordinary least-squares, compressive sensing and least-angle- 
regression. 

Least angle regression (LAR) ETll is an iterative regression method for variable selection, 
which was proven to be particularly powerful when the number of regressors is much larger 
than the number of experiments, and only a few of them are relevant. Assuming that the model 
of interest respects the principle of sparsity, i.e. it can be approximated using a relatively small 
number of polynomials, LAR is used to derive sparse polynomial chaos expansions (SPCE) in 
H. At each iteration, a basis set is retained with the associated coefficients and a measure of 
accuracy, e.g. the leave-one-out error, is computed. At the end of the process, based on the 
accuracy measure, one can choose the best performing PCE model. Note that after detecting 
the relevant basis by LAR, one should recompute the coefficients using ordinary least-square 
method (hybrid-LAR) E4ll . The LAR algorithm can be summarized as follows: 

1. Initialize the set of candidate regressors to the full basis and the set of selected regressors 
to 0. 
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2. Initialize all coefficients equal to 0. Set the residual equal to the output vector. 

3. For each iteration until the stop condition is satisfied or all the regressors have been ana¬ 
lyzed, perform the following steps: 

• Determine the most correlated candidate basis element to the current residual and 
add it to the selected basis 

• Move simultaneously all the coefficients associated with the selected basis until 
another basis element is as correlated with the residual as they are. Update the 
residual. 

4. Choose the best iteration based on a prescribed error estimate, e.g. eLoo- 

Readers are referred to for further details. The algorithm is schematically illustrated in 
Figure [TJ in which the adaptive selection of the relevant basis is presented for the first 5 terms 
in the expansion. The SPCE based on low-rank hyperbolic truncation and LAR is hereafter 
considered as the reference approach. 
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(b) Iteration 1 
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(f) Iteration 5 


Figure 1: Least angle regression in the multi-index space, p = 5, r = 2 and g = 0.75 are used for the truncation 
schemes. The cross symbols (x) represent the candidate basis. The dot symbols (•) represent the basis selected by 
LAR in the current iteration. 
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3 HIERARCHICAL ADAPTIVE POLYNOMIAL CHAOS EXPANSIONS 

Recently, ifTOll and ifTTll proved that adaptive basis scheme can be more efficient than fixing 
a candidate basis a priori. Moreover, ifTOll propose a novel adaptive-basis non-intrusive spec¬ 
tral projection-based approach for PCE. An initial basis is enriched by adding candidate sets 
associated with adaptive sparse grids and removing the unimportant PC basis elements with 
minor contributions to the variance. In IfTTll . it is assumed that a basis function is admissible, or 
potentially relevant, only if its backwards neighbors exist in every dimension. Given the same 
computational budget, the authors observe that adaptive basis selection schemes provide more 
accurate PCE than predefined basis schemes. 

In this section, we introduce a novel adaptive basis SPCE scheme which relies on LAR and 
the so-called principle of heredity. We emphasize that the models of interest are also assumed 
to follow the principle of sparsity. 

3.1 Principle of heredity 

Heredity is a biological process in which the characteristics of parents are passed down to 
their child through the genes. This process is represented by Mendel’s principle of heredity 
ifTSll which was recently generalized in a statistical sense in lfT9ll . According to the principle 
of heredity, if two factors are not relevant in the model, it is likely that their interaction is also 
irrelevant. This principle is nowadays used for the detection of important interaction effects 
(interaction screening) in high and ultra-high dimensional problems. It was incorporated in 
variable selection techniques, e.g. forward selection scheme ll 20 ll . LAR ifTSll and structured 
variable selection and estimation lIlTIl . 

Herein, we introduce the use of the heredity principle in the context of SPCE. Consider a 
the model Y{X) with X denoting the vector of M independent input parameters. Based on 
the hierarchy principle, one assumes that Y can be approximated with only main effects and 
two-dimensional interaction effects as follows: 

Y = yo+ ypxpp{X) + y~^'ip~^{X) + e (8) 

p&Ai 76.42 

where Ai = {f3 : ||/3||o = 1, ||/3||i ^ p} and A = {7 G : ||7||o = 2 , || 7 ||i ^ p}, 

i.e. the second (resp. third) summand consists of univariate (resp. bivariate) polynomials; yo is 
the constant term, y/^ and y~^ are the coefficients of the expansion and e is the truncation residual 
error. 

The heredity principle states that the interaction term ip~^{X) = 7 ^ j is 

active in the model only if the parent terms ff{Xi) and fi^.{Xj) are also present EOll . More 
precisely, two forms of the heredity principle may be considered, namely weak heredity and 
strong heredity: 

• Weak heredity indicates that given an active interaction term {Xi)fi^. (Xj), at least one 
of the two parent terms (Xi) and (Xj) exist in the model. 

• Strong heredity implies that both parent terms (Xi) and (Xj) must be active for the 
interaction term {Xi)ipi^. (Xj) to be included in the model. 

Figure illustrates the principle in the example of a four-dimensional model, which requires 
only first degree univariate polynomials and their interactions. 
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(a) Weak heredity 


(b) Strong heredity 


Figure 2: Heredity principle 


To introduce the principle of heredity in SPCE, we need to combine this principle with the 
LAR method. To this end, ifTSll introduced distinct generalized LAR algorithms for strong and 
weak heredity. The proposed algorithms, however, cannot solve problems in which both strong 
and weak heredity are active simultaneously, e.g. one group of variables follows strong heredity, 
another group follows weak heredity. A clear distinction between strong and weak heredity in 
a practical problem is, however, quasi-impossible a priori. This limitation is overcome by the 
approach proposed in the next section. 

3.2 Hierarchical adaptive PCE 

An iterative approach, named hierarchical adaptive PCE, is proposed to incorporate both 
forms of heredity in the LAR scheme. After detecting the most relevant term in the current 
iteration, the heredity principle is used for updating the current candidate set (which is used for 
later iterations). The proposed scheme works as follows: 

1. Generate an initial candidate set, which consists of one-dimensional (1-D) polynomial 
chaos basis 

2. Eor each iteration, perform the following steps until the stop condition is satisfied: 

(a) Find the most relevant term in the current candidate set. 

(b) If it is a 1-D term 

i. Generate interaction terms between the 1-D term selected in this iteration and 
the remaining 1-D terms in the current candidate set. 

ii. Compare the relevance of the selected 1-D term with the generated interaction 
terms. The relevance of a basis element is represented by the correlation coef¬ 
ficient between it and the current residual. If the selected 1-D term remains the 
most relevant, then: 

• Compute the associated coefficients as in LAR 

• Update the candidate set by adding the interaction terms generated in the 
current iteration. 

iii. otherwise, select the most relevant interaction term, compute the associated 
coefficients as in LAR 
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(c) otherwise, it is an interaction term. Compute the associated coefficients as in LAR. 

(d) Extract the selected term from the candidate set. Update the residual. Check the 
stop condition. 

The approach is illustrated in Figure until six terms are active in the expansion. In iteration 4, 
the interaction term (ai, 02 ) = (2, 3) is selected after its parent terms (2,0) and (0, 3) have been 
selected in the previous iterations. Thus, the selection in the current iteration follows the strong 
heredity principle. In iteration 5, the I-D term (5, 0) is selected first, but it is less relevant than 
its child interaction term (5,3), which is generated on-the-fly. The latter is therefore selected 
instead of its parent term. Only one parent of this interaction term, i.e. (oi, 02 ) = (0,3), is 
currently active. This selection illustrates the weak heredity principle. In iteration 6, the weak 
heredity prevails again when the interaction term (2,4) is selected instead of its parent (0,4). 

Now we explain the superiority of the proposed approach compared to ifTSll . in which authors 
proposed distinct algorithms for strong and weak heredity. If only the strong heredity principle 
is used, then 1-D term s must always be selected before their interaction terms. This may intro¬ 
duce error, when the interaction effects are important ll22l . If one uses only the weak heredity 
(which already covers the strong heredity), the interaction terms would quickly include the full 
candidate basis. In addition, in practical problems, it is impossible to know in advance whether 
the considered model follows the heredity principle or not, and if it does, in which form. In 
particular, a group of variables may follow the strong form whereas another group the weak. 
Our approach is flexible in the sense that both forms of heredity principle can be considered 
simultaneously in each iteration. 

It is worth to emphasize that one can also apply the hyperbolic and low-rank truncation 
schemes when generating the interaction terms according to the heredity principle. This helps 
reducing further the size of the candidate set. In this case, the proposed approach can be consid¬ 
ered as an extension of the reference approach except that the candidate polynomial chaos basis 
is updated iteratively by means of the principle of heredity. The size of the candidate set in the 
proposed algorithm increases iteratively, while remaining always smaller than that of the refer¬ 
ence approach. Therefore, it is expected to outperform the the reference approach, in particular 
when the considered model requires anisotropy in the polynomial basis, i.e. high polynomial 
degrees in certain directions versus low degrees in others. 

As a summary, this approach allows one to update the candidate basis in each iteration using 
the heredity principle rather than using a fix candidate set. One can discard a large number of 
interaction terms that are not relevant and just focus on the significant ones that are selected 
on-the-fly. 

4 NUMERICAL EXAMPLES 

The proposed hierarchical adaptive PCE (hereafter denoted by h-LAR) is implemented in 
UQLab ll2^ (the Uncertainty Quantification toolbox in Matlab) developed at the Chair of Risk, 
Safety & Uncertainty Quantification (ETH Zurich, Switzerland), which already includes the 
reference approach (hereafter denoted by LAR). 

Let us examine the effectiveness of h-LAR with respect to LAR by means of two numerical 
examples. 
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Figure 3: Hierarchical adaptive PCE based on the heredity principle in the multi-index space, p = b and r = 2 are 
used for the truncation scheme. Second-rank PCE are assumed to be sufficiently accurate. The cross symbols (x) 
represent the candidate basis. The dot symbols (•) represent the basis selected in the current step. 


4.1 SoboP function 

We consider the Sobol’ function llTdll : 


y=Il 


|4Xj — 21 + Q 

1 + Ci 


(9) 


in which c = {1, 2,5,10, 20, 50,100,500}^ and the input parameters are uniform random vari¬ 
ables Xj ~ ^Y[0,1], i = 1, ... , 8. The Sobol’ function is highly anisotropic in that it requires 
high polynomial degree for low-i dimensions. 
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We first investigate the effect of the experimental design’s size on the accuracy of the result¬ 
ing PCE models. For a given size (e.g. 100,150, 200, 300,400, 500,1000), one repeats the Latin 
hypercube sampling (LHS) of the experimental design 100 times then builds the correspond¬ 
ing PCE by means of LAR and h-LAR. The mean value and 95% confidence interval of the 
estimated LOO errors are depicted in Figure]^ Note that LOO error is chosen as an indicator 
of accuracy because it is much more sensitive to overfitting than the commonly used empiri¬ 
cal error m. h -LAR performs significantly better than LAR with faster convergence rate and 
smaller variability due to LHS. 

As a further validation, given an ED of size 200, we compute the PCE by the two approaches 
using the same truncation options (p = 30, g = 0.5, r = 2). Using the computed PCE, we 
predict the output values on a new set of input parameters {A'j, i = 1, ... , 10^}. The predicted 
values of the output and its probability density function (PDF) are compared with the values 
and the PDF computed with the analytical model. The results are shown in Figure]^ Similarly, 
Figure depicts the results when an ED of size 300 is used. The corresponding results (e.g. 
polynomial degree giving the best accuracy, relative validation error ey and number of retained 
polynomials Nj.) of the PCE are given in Table [Tj In both cases, h-LAR outperforms LAR in 
predicting specific values of the output as well as its statistical distribution. One sees that given 
the same computational budget (i.e. the same ED), h-LAR allows one to reduce the validation 
error by a factor of 5. The best degree in h-LAR is higher than that in LAR, which leads to a 
larger number of selected terms in the expansion. 



Figure 4: Mean value and 95% confidence interval of the leave-one-out error as a function of the experimental 
design’s size. The mean values and confidence intervals are obtained with 100 Monte Carlo realizations. 


Method 

Size of ED 

Q 

r 

Best degree 


Nr 

LAR 

200 

0.5 

2 

7 

5.95 X 10-2 

37 

h-LAR 

200 

0.5 

2 

9 

1.6 X 10-2 

100 

LAR 

300 

0.5 

2 

9 

2.2 X 10-2 

88 

h-LAR 

300 

0.5 

2 

20 

4.3 X 10-3 

154 


Table 1: Options and results of the utilized PCE 
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(a) Low-rank hyperbolic truncation scheme 



(b) Hierarchical adaptive PCE 

Figure 5: Sobol’ function - Experimental design of size 200 - PDF of the output. 


4.2 Schwefel function 

We now consider a modified version of the Schwefel function ll^ which reads: 


(^ + 0.5)X,sin(^^ |X,|) 

p- ( 10 ) 

+ ^—V(—+ 0.5)W sin(w2- |X,|) VXfc 

3000 ^^20 ^ ^ ^V20' 

jG*Si k^S2 

where Xi - f/[-500,500], i = 1, ... ,20, Si = {1,3,5,6,8} and .Sa = {15,18,20}. This is 
a complex function with many local minima. The first term is modified in such a way that the 
function requires high polynomial degree in high-z dimensions and the second term is added to 
the function in order to introduce interaction effects between some of the dimensions. 

Given an ED of size 1,000, we compute the PCE of the Schwefel function by means of the 
two approaches with the same truncation options, i.e. p = 30, q = 0.25, r = 2. The values 
of the function on 10^ samples of the input parameters are predicted by the two PCE models 


i=l 
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(a) Low-rank hyperbolic truncation scheme 



(b) Hierarchical adaptive PCE 

Figure 6: Sobol’ function - Experimental design of size 300 - PDF of the output. 


and compared with the values obtained with the numerical functions. The outcomes of the PCE 
are presented in Table The relative validation error of h-LAR is smaller than that of EAR, 
whereas more polynomials are retained by h-LAR. For this numerical model, h-LAR performs 
slightly better than LAR. 


Method 

Size of ED 

q 

r 

Best degree 


Nr 

LAR 

1000 

0.25 

2 

23 

2.2 X 10-2 

556 

h-LAR 

1000 

0.25 

2 

30 

1.03 X 10-2 

671 


Table 2: Options and results of the utilized PCE 


5 CONCLUSIONS AND PERSPECTIVES 

In modem engineering and computational sciences, polynomial chaos expansions (PCE) 
have been widely used as a powerful tool for uncertainty quantification. The application of 
PCE to high-dimensional complex problems might be, however, hindered due to computational 
limitations. The reference sparse PCE technique relies on selecting the relevant polynomial ba- 
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sis functions from a candidate set defined a priori by means of appropriate truncation schemes. 
Recent research shows that it can be more effective to compute PCE with an adaptive basis 
enriehment approaeh instead. 

In the current paper, we propose a novel hierarehical adaptive-basis approach, which con¬ 
sists in selecting the important basis elements from an iteratively enriched candidate basis. The 
proposed approach is based on combining the least-angle-regression method with the principle 
of heredity. The latter, which originated in a biological context, is herein utilized in a statis¬ 
tical sense for detecting the potentially relevant interaction effects in the model. Using some 
numerieal models, we proved the effectiveness of our approaeh. 

Further investigations are required in order to take into account higher interaction order. In 
addition, the current approach updates the candidate basis set only by adding new terms, while 
it can also be important to remove irrelevant terms during the process, which may contribute to 
further deereasing the computational costs and increasing the aeeuraey of PCE. 

6 Acknowledgment 

The first author would like to thank Dr. Stefano Marelli (ETH Ziirieh, Chair of Risk, Safety & 
Uncertainty Quantifieation) for various fruitful discussions during the preparation of the paper. 
We also acknowledge the Development team of UQLab (the Uncertainty Quantification toolbox 
in Matlab) for providing the numerical platform on which this work is carried out. 


REFERENCES 

[1] D Xiu and G E Kamiadakis. Modeling uncertainty in flow simulations via generalized 
polynomial chaos. J. Comput. Phys., 187:137-167, 2003. 

[2] M Paffrath and U Wever. Adapted polynomial chaos expansion for failure detection. J. 
Comput. Phys, 226(1):263-281, 2007. 

[3] K Young and M Grace. Simulation of stochastic quantum systems using polynomial ehaos 
expansions. Phys. Rev. Lett., 110(11): 110402, 2013. 

[4] R H Vaz, J M C Pereira, A R Ervilha, and J C F Pereira. Simulation and uneertainty 
quantification in high temperature microwave heating. Appl. Therm. Eng., 70(1): 1025- 
1039, 2014. 

[5] S K Choi, R V Grandhi, R A Canfield, and C L Pettit. Polynomial ehaos expansion with 
Latin Hypercube sampling for estimating response variability. AIAA Journal, 45:1191- 
1198,2004. 

[6] M Berveiller, B Sudret, and M Lemaire. Stochastic finite elements: a non intrusive ap¬ 
proach by regression. Eur. J. Comput. Mech., 15(l-3):81-92, 2006. 

[7] A Doostan and H Owhadi. A non-adapted sparse approximation of PDEs with stoehastic 
inputs. J. Comput. Phys, 230(8):3015-3034, 2011. 

[8] L Fagiano and M Khammash. Nonlinear stoehastic model predietive control via regu¬ 
larized polynomial chaos expansions. IEEE 51st Annual Conference on Decision and 
Control (CDC), pages 142-147, 2012. 


13 




Chu V. Mai, Bruno Sudret 


[9] G Blatman and B Sudret. Adaptive sparse polynomial chaos expansion based on Least 
Angle Regression. J. Comput. Phys, 230:2345-2367, 2011. 

[10] Z Perko, L Gilli, D Lathouwers, and JL Kloosterman. Grid and basis adaptive polynomial 
chaos techniques for sensitivity and uncertainty analysis. J. Comp. Phys., 260:54-84, 
2014. 

[11] J D Jakeman, M S Eldred, and K Sargsyan. Enhancing /i-minimization estimates of 
polynomial chaos expansions using basis selection. ArXiv:1407.8093, 2014. 

[12] C Soize and R Ghanem. Physical systems with random uncertainties: chaos representa¬ 
tions with arbitrary probability measure. SIAM J. Sci. Comput., 26(2):395^10, 2004. 

[13] D Xiu and G E Kamiadakis. The Wiener-Askey polynomial chaos for stochastic differen¬ 
tial equations. SIAM J. Sci. Comput., 24(2):619-644, 2002. 

[14] G Blatman. Adaptive sparse polynomial chaos expansions for uncertainty propagation 
and sensitivity analysis. PhD thesis, Universite Blaise Pascal, Clermont-Eerrand, 2009. 

[15] M Yuan, V R Joseph, and Y Lin. An efficient variable selection approach for analyzing 
designed experiments. Technometrics, 49(4):430^39, 2007. 

[16] B Sudret. Polynomial chaos expansions and stochastic finite element methods. In K K 
Phoon and J Ching, editors. Risk and Reliability in Geotechnical Engineering, chapter 5. 
Taylor and Erancis, 2014. 

[17] B Efron, T Hastie, I Johnstone, and R Tibshirani. Least angle regression. Annals of 
Statistics, 32:407^99, 2004. 

[18] W Bateson and G Mendel. Mendel’s principles of heredity. Courier Corporation, 2013. 

[19] M Hamada and C E J Wu. Analysis of designed experiments with complex aliasing. J. 
Qual. TechnoL, 24(3):130-137, 1992. 

[20] N Hao and H H Zhang. Interaction screening for ultra-high dimensional data. J. Amer. 
Stat. Assoc., 109(507): 1285-1301, 2014. 

[21] M Yuan, V R Joseph, and H Zou. Structured variable selection and estimation. The Annals 
of Applied Statistics, pages 1738-1757, 2009. 

[22] R D Shah. Modelling interactions in high-dimensional data with backtracking. Technical 
report. Statistical Laboratory, University of Cambridge, 2012. 

[23] S Marelli and B Sudret. UQLab: a framework for uncertainty quantification in Mat- 
lab. In ASCE, editor, 2nd Int. Conf. on Vulnerability, Risk Analysis and Management 
(ICVRAM2014), pages 2554-2563, 2014. 

[24] M Laguna and R Marti. Experimental testing of advanced scatter search designs for global 
optimization of multimodal functions. J. Global Optim., 33(2):235-255, 2005. 


14 




